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ABSTRACT 


Using the geometrical theory of diffraction, the traditional 
method of high frequency scattering analysis, the prediction of 
the radar cross section of a perfectly conducting, flat, 
rectangular plate is limited to principal planes. Part A of this 
report predicts the radar cross section in nonprincipal planes 
using the method of equivalent currents. This technique is based 
upon an asymptotic end-point reduction of the surface radiation 
integrals for an infinite wedge and enables nonprincipal plane 
prediction. The predicted radar cross sections for both 
horizontal and vertical polarizations are compared to moment 
method results and experimental data from Arizona State 
University's anechoic chamber. 

In part B of this report, a variational calculus approach to 
the pattern control of the horn antenna is outlined. The approach 
starts with the optimization of the aperture field distribution so 
that the control of the radiation pattern in a range of directions 
can be realized. A control functional is thus formulated. Next, 


a spectral analysis 

method is 

introduced 

to 

solve for the 

eigenfunctions from 

the extremal 

condition 

of 

the formulated 


functional. Solutions to the optimized aperture field 


distribution are then obtained. 


Fart A 


Nonprincipal Plane Scattering of Flat Plates 
I. Introduction 

The design of low-observable vehicles, with a reduced radar 
cross section (RCS), is a problem of current interest within the 
scientific and military communities. The parameters of interest 
in RCS reduction are the material composition and shape of the 
target. The importance of each of these factors varies with the 
type of radar used for detection. Over-the-horizon radar, used 
for long-range surveillance, operates at low frequencies in the 
VHF region or lower, resulting in targets on the order of one or 
two wavelengths in size. The material composition of the target 
is the chief factor in RCS minimization in this frequency range. 
Other radars, such as airborne systems, operate at higher 
frequencies, usually above 1 GHz, rendering targets several 
wavelengths in size. In this frequency range the shape of the 
target takes precedence over material composition in RCS 
determination. 

For high-frequency scattering prediction, each part of the 
body of interest can be viewed as an independent scatterer. RCS 
prediction at high frequencies thus involves decomposing the 
target into simple geometric entities, calculating the RCS of each 
entity, and then appropriately summing the individual predictions 
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to achieve an overall RCS. The accurate modeling of simple shapes 
is, therefore, critical. One geometry of particular interest is 
the perfectly conducting, flat, rectangular plate. 

The Geometrical Theory of Diffraction (GTD) fl] and the 
Uniform Theory of Diffraction (UTD) [2], the traditional 
high-frequency modeling techniques, accurately predict the 
scattering from the flat plate in the principal plane away from 
grazing and near-grazing incidences [3]. The GTD/UTD results are 
inaccurate in directions away from the Keller cone of diffracted 
rays; therefore, this method cannot be used in modeling the 
scattering patterns in nonprincipal planes, or planes not parallel 
to one of the edges of the plate. Ross's analysis [3] of the 
principal plane scattering of the plate involves truncating the 
two-dimensional GTD solution for an infinite, perfectly 
conducting, finite-width strip. This truncation does not account 
for the presence of the two edges parallel to the plane of 
incidence and diffraction. These edges are not critical in 
scattering near normal incidence; but as the grazing and 
near-grazing directions are approached, the role of these two 
edges in the scattering mechanism becomes more important. Because 
the GTD solution does not include these edges, the predictions 
near grazing and at grazing do not agree with experimental 
results . 

Sikta [4] developed another high-frequency analysis using the 
method of equivalent currents. By dividing the plate into strips 
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that approach zero thickness as a limit and summing the 
contributions from each strip, the RCS for principal and 
nonprincipal scattering planes can be calculated by varying the 
orientation of the strips. The RCS of each strip is calculated 
using the method of equivalent currents. 

The method of equivalent currents lacks the simplicity of the 
GTD approach to scattering prediction and the insight into 
physical scattering points that GTD provides. Despite these 
drawbacks, it has two distinct advantages over GTD. First, it can 
be used to correct for axial caustics encountered by GTD [5]. 
Second, and more importantly for the modeling of the flat 
rectangular plate, the method of equivalent currents can be used 
to predict scatttering in directions away from the Keller cone. 

The general theory behind the method of equivalent currents 
is that the fields scattered by an object are approximated by the 
fields radiated by nonphysical electric and magnetic currents 
placed along the edges of the scatterer. The crucial problem in 
scattering prediction using the method of equivalent currents is, 
therefore, the formulation of accurate equivalent currents. 
Several different versions of equivalent currents exist. One of 
the first references to the method of equivalent currents is in 
Millar's work on scattering by a circular aperture in a perfectly 
conducting, plane, infinite screen with plane wave incidence [6]. 
These equivalent currents are derived by extracting the first term 
in the asymptotic expansion of the exact Sommerfeld half-plane 
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solution. The extent of the applications of Millar’s equivalent 


currents is limited. 

Ryan and Peters formulated equivalent currents using the GTD 
diffraction coefficient [5]. Conceptually these are the same as 
Millar's currents, although they differ in form. These 
successfully correct for axial caustics [5] and predict radiation 
patterns of rectangular waveguides [7] and horn antennas [ 8 J . 
Because of the nature of the GTD diffraction coefficient, their 
equivalent currents are limited to scattering prediction on the 
Keller cone. Despite this limitation, Sikta predicts nonprincipal 
plane scattering from flat plates using these equivalent currents 
by careful orientation of the strips that divide the plate. 

Michaeli formulated more rigorous equivalent currents that 
are valid for scattering directions away from the Keller cone of 
diffracted rays [9-11]. These currents are derived using an 
asymptotic end-point reduction of the surface radiation integral 
to a line radiation integral. The GTD equivalent currents [9] 
contain infinite discontinuities at certain directions of 
incidence and observation. By considering the integration across 
the surface at a skew angle, the Physical Optics/Physical Theory 
of Diffraction (PO/PTD) equivalent currents [10-11] remove all 
singularities, except the Ufimtsev singularity for forward 
scattering at grazing incidence. 

This work examines the monostatic RCS of the perfectly 
conducting, flat, rectangular plate using the GTD [9] and PO/PTD 



[10-11] equivalent currents of Michaeli. Both horizontal and 
vertical polarizations are considered. Scattering in nonprincipal 
planes is emphasized. Results are compared with data from moment 
method (MM) computations and experimental measurements performed 
in Arizona State University's (ASU's) compact range. 



1 1 . Theory 


A. RCS Calculation 


The plate geometry, shown in Fig. J, is oriented so that the 

x-z and y-z planes correspond to principal planes. Nonprincipal 

planes are designated by a tilt angle 4 >’ , measured from the 

positive x axis. RCS values in a designated plane are calculated 

jot 

as a function of 0' , measured from the positive z axis. An e 
time convention is assumed and suppressed throughout. 

The standard definition of the three-dimensional RCS is 


[ 12 ]: 


o 


3-D 


=lim 

r->® 




} 


( 1 ) 


This definition implies that the observation point is in the 

2 

far-field region, or a distance of at least 2D A from the target 
where D is the largest dimension of the target. 

Both horizontal and vertical polarizations are considered. 
For the horizontal polarization the incident electric and magnetic 
fields are: 


H.= a 0 E o (l/n)e‘ j -°- 

The fields for vertical polarization are: 

A • i k ° r 

V a e V - - 

V Vo (1/ri)e ' j " 


(2a) 

(2b) 


(2c) 

(2d) 
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a^ and a~ are the standard directional unit vectors in the 

<P 0 

spherical coordinate system, and k is the wavenumber, 2 n/X, in the 
direction of propagation. 

In order to simplify the analysis for the rectangular plate, 
the incident fields will be transformed to the rectangular 
coordinate system. The position vector £ is: 


r~ax + ay + az 
— x y z 


(3) 


The propagation vector k for the incident field is: 


k = -ka = -k(a sin0cos0 + a sin0sin0 + a cos0) (4) 

— r x y z 

With respect to the rectangular coordinate system, the incident 
fields are: 

Horizontal Polarization 


E. =E e 

1 o 


j k ( x s i n0co 8 0+ y e i n0e i n<J>+ zc o 8 0 ) 


[-a sin0 +■ a cos<J> ] (5a) 

x y 


, „ . v j k ( x s i nOco o0+y s i n08 i n0+zco s0 ) w 

H.=E (l/rj)e X 

1 o 

[a cos0cos0 + a cos0sin0 + a sin0 ] 

x y z 


Vertical Polarization 


j k ( xe i n0cos0+y8 i n0s i n0+zcoc0) 

E — E e X 

i o 

[a cos0cos0 + a cos0sin0 + a sin0 ] 
x y z 

„ T , . „ , . j k ( xs i n0cos0+ys i n0s i n0+zco80) 

H — E (1 /n ) e X 

~ i o 

[-a sin0 + a cos<J> ] 

x y 


(5b) 


(5c) 


(5d) 


To calculate the scattered fields, a vector potential 
approach is used. The electric field components in the far zone 
are approximated by [13]: 


E ss 0 

r 


(6a) 
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V -» 

1*6 + "V 

(6b) 

V ~ j “ 

w. - 1F 0 ] 

(6c) 

The magnetic (A) 

and the electric (F) 

vector potentials are 

defined as [ 13] : 



J w 

ii 

<1 

I(x’ , y ' , z ' ) ^ dl * 

(7a) 

F - ~ 

- 471 J 

c 

e‘ jkK 

M(x ' ,y 1 ,z' ) ^ dl ' 

(7b) 


X and M are the equivalent electric and magnetic currents, 
respectively, which are placed along the perimeter of the 
scatterer. They radiate into free space according to (6) and (7) 
and approximate the scattered fields. The geometry for far-field 
scattering is shown in Fig. 2. The location of the scatterer is 
designated by primed coordinates and the observation location by 
unprimed coordinates. Integration is along C, the edge of the 
scatterer. The distance from the source point to the observation 
point is R, represented by 

R=|R|=|r-r'| (8) 

The location vectors of the source and observation points with 
respect to the origin of the coordinate system are r/ and jr, 
respectively . 

For far-field scattering the following simplifications can be 

made : 

R £ r-r’cos# = r - a (for phase variations) (9a) 

R £* r (for amplitude variations) (9b) 

For the flat plate centered at the origin of the rectangular 


coordinate system: 



r ' = a x + a y 

— x y 


( 10 ) 


Thus, for phase variations: 

R as r-xsinOcos<f> - ysinOsin<f> (11) 

By substituting (9b) and (11) into (7a) and (7b), the vector 


potentials for the flat plate become 
•jkr f 


v= iL iL 


- 4n r 

- jkr 


F= — — 
- 471 r 


_ / . j k ( x b i n0co 80+ y 8 i n0e i r»0) 

_I(x,y,z)e J dl 


XMt v j k ( xs i n0co80+y e i n08 i n0) , ^ 

M(x,y,z)e dl 


(12a) 


(12b) 


These integrals are evaluated along the perimeter of the plate. 
To simplify the derivation, each integral is represented as a sum 
of four integrals, each corresponding to an edge of the plate. 
For the plate the integration is of simple form, and the 
potentials can be expressed in closed form. 1^ anc * M must be 
determined separately for each edge taking into account the 
individual geometries. Each edge is viewed as the truncation of 
an infinite wedge, which reduces to a half plane for the flat 
plate. The formulation for a general half plane follows. 


B. Equivalent Currents 


1. Geometry 

The general wedge configuration used for both the GTD [9] and 
PO/PTD [10-11] equivalent currents is shown in Fig. 3. For a half 
plane, n, which denotes the wedge angle, is 2. There are several 
directional vectors and angles which are crucial to the 
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development. The directional vectors are: 

s 1 — unit vector in the direction of Incidence 

s= unit vector in the direction of observation 

t= unit vector tangent to the edge of interest, directed 
so that it encircles the scatterer in a 
counterclockwise manner 

n= unit vector normal to the edge of interest, lying on 
the upper face 

The angles are: 

p'= angle between s' and the edge 
p = angle between s and the edge 

'9'= angle between the upper face and the edge-fixed 
plane of incidence 

d = angle between the upper face and the edge-fixed 
plane of observation 

a = skew angle for integration across the surface 


In terms of the directional vectors, the angles are: 

P ' = cos 1 (s ' o t ) (13a) 

p = cos’ 1 (s<>t) (13b) 


•9 ' = cos 


-l 


s'x t 
|s'X t| 


X t 


o n 


■9 = cos 


-1 


s x 4- x ; 


|s X t| 


:} 


(13c) 


(13d) 


This report considers only backscattering from the flat, 
rectangular plate. For backscattering, 0=Jr-/3’ , and s=~s 1 . 

For the flat plate this vector is: 

s’= -a = -a sin0cos4> - a sin0sin<#> - a cos0 (14) 

r x y z 

A different set of angles must be formulated for each edge. 
The vectors for each edge are: 
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Edge 1 : 


t = a 
i y 


n = -a 

1 X 


Edge 2: t = -a 

2 y 


Edge 3: t = -a 

3 x 


n - a 

2 x 


n = -a 
3 y 


1 = a 

A x 


n - a 

A y 


Edge 4: 

For each edge the tangential components of the incident 
electric and magnetic fields are needed to determine the 
corresponding equivalent currents. The tangential electric field 


component at the edge of interest is: 

c ian c I 
E * E ° t 
i i 

and the magnetic field component is: 
H tan = H o t 

i “i 


(15a) 


(15b) 


The individual current components are easily determined by 
substitution of the incident field components and directional 
angles into the appropriate current equation. 


2. GTD Equivalent Currents 

The GTD equivalent currents [9] are based upon an asymptotic 
endpoint evaluation of the surface radiation integrals for a 
wedge. The integration across the surface is along a coordinate 
normal to the edge of the wedge. This choice of integration 
coordinate results in infinite discontinuities in the calculated 
fields when bistatic calculations are made. The monostatic 
results using these currents, however, avoid infinities. 

The expressions for the equivalent electric (I) and magnetic 



(M) currents for a general wedge configuration are quite long and 
are, therefore, not repeated here. For backscattering 
fi=n -0, and s=-s') from a half plane (n=2), they simplify 
considerably, and they take the forms of: 

Electric Equivalent Current : 

1 J tan (cosfl-1 ) _ „tan 2COS0 1 ( l+cosfl ) 

.. 2 0l | i r/cosO i sind 

jksin 0 ' t 

Magnetic Equivalent Current : 

M= H tan nill£°s*± (16b) 

1 jksin 2 p'cosd 

An interesting aspect of these currents is that the electric 
current depends on both the incident tangential electric and 
magnetic fields. Earlier versions formulated by Millar [6] and 
Ryan and Peters [5] do not contain this dependence, which makes 
them invalid in directions away from the Keller cone of diffracted 
rays. The currents of (16a) and (16b) are valid in all directions 
of observation. 

3. PO/PTD Equivalent Currents 

In order to correct for the discontinuities incurred by the 
GTD equivalent currents, Michaeli derived new currents [10-11] 
that incorporated integration across the wedge at a skew angle 
rather than in a direction normal to the edge under consideration. 
The skew coordinate of integration corresponds to the grazing 
diffracted ray, which is shown in Fig. 3. The angle of this ray 


(16a) 
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with respect to the edge is p' t the same as the angle that the 
incident ray makes with the edge. This choice of geometry for the 
integration eliminates all singularities except the Ufimtsev 
forward scatter singularity due to grazing incidence. 

The PO/PTD analysis consists of PO and PTL) current 
components. The PO components [10] account for the radiation of 
the surface currents that would be present if the scatterer were 
infinite in extent. The PTD components [11] account for the 
finite boundaries of the scatterer. The sum of these components 
is finite in all directions except for the Ufimtsev singularity 
for forward scatter and grazing incidence. 

A. PO Components 

The forms of the general PO equivalent currents [10] are 
rather tedious expressions involving the modified Fresnel 
transition function. For far-field, plane-wave backscattering 
from a half plane, however, the Fresnel transition function can be 
replaced by its asymptotic expansion. In this case the components 
simplify to: 

po J sinasinfl* ~ _ * 1 

kr] sin£' — i (sin0sin£ ’ costf * - cos0cos£’) & 


„po 5 _ _ siriasin# * ~ ) 1 

k sin#’ 1 — i ' (sin0sin£ ' cos# ' - cosOcosjS 1 ) 


(17b) 


The angle a denotes the direction of the transverse 
coordinate for integration. For the analysis in this work, a was 

o 

taken to be 90 , corresponding to integration normal to the plate 
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o 

edge. Accurate rc?sults are obtained for a between <P and 90 . 
B. PTI) Components 

The PTD components that must be added to the PO currents 
reduce to the following for backscattering from a half plane: 


r P t d 


,tan 


sin(9 1 /2 ) 


j/2 

ksin^ ^ * r)sin0 ' (cos9 ' -cot^f} ' ) 

[]/ l-cos9 ’ +2cot 2 0 ' /2 cos(9'/2) ] + 

j|tan 2cos ( 9 1 /2 ) cotp ' 

J l-cos9 ' +2cot Z 0 ' 
jr}sin9 ' 


M p td m H Un 


! ksin Z j5’ (cos9 '-cot 2 0' ) 
/2 cos (9 ' /2 ) 


1-cos 


cos-9 ' +2cot 0 


(18a) 


(18b) 


Note that, similar to the GTD electric currents, the PTD electric 
currents depend on both the incident electric and magnetic fields. 


C. Assembly of Scattered Fields 

The calculation of the RCS involves the assembling of the 
appropriate equations. First, calculate the currents for each 
edge of the plate using either the GTD currents of (16a) and (16b) 
or the PO/PTD currents of (17a)~(18b). To determine the 

tangential fields, use (15a) and (15b) and either (5a) and (5b) 
for horizontal polarization or (5c) and (5d) for vertical 
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polarization. Substitute the appropriate currents into (12a) and 
(12b) to obtain the vector potentials. Finally, substitute into 
(6b) and (6c) to obtain the electric field components. 

Unlike the fields predicted by traditional GTD analysis, the 
scattered electric fields for both polarizations (as predicted by 
the proposed methods) include cross polarization terms. However, 
in the calculation of the RCS, only the primary polarization terms 
are considered. In other words, the components are used for 
horizontal polarization, and the Eg components are used for 
vertical polarization. 

A final note of interest in the RCS calculations using the 
equivalent currents method is that the fields must be multiplied 
by a factor 1/2 to obtain the correct magnitudes. This is a 
problem that Ryan and Peters [5] and Sikta [4] also encounter. 
Ryan and Peters were off by a factor of 1/2 elsewhere in their 
derivation and, therefore, did not recognize that the fields were 
double the correct values. Sikta pointed out their error and 
avoided the same problem by considering currents only along edges 
normal to the plane of incidence and diffraction in calculating 
the fields. One intuitive explanation for the necessity of the 
factor of 1/2 is that in placing currents along all four edges of 
the plate, the integration across the plate is performed twice, 
thus doubling the fields. For a four-sided parallelogram the 
compensating factor of 1/2 is easily justified. It is not clear 
what the correction factor should be in applying the equivalent 



currents method to scatterers with more edges. This will be 
examined in our next period. 
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III. Results 


Computations were made for a square plate with each side 
equal to 5.73 X. The results of both the GTD and PO/PTD 
equivalent currents analyses were compared with MM computations 
and measurements performed at ASU's anechoic chamber. The results 

o 

for scattering in a plane designated by <f>'=30 are shown in Fig. 4 
for horizontal polarization and in Fig. 5 for vertical 
polarization. Figs. 6 and 7 are for horizontal and vertical 

o 

polarization, respectively, at a tilt angle of 4>'=45 . 

The GTD equivalent currents results agree very well with the 
PO/PTD results, which is expected for monostatic calculations. 
For bistatic computations, the GTD currents encounter infinities 
in many directions. The PO/PTD formulation avoids infinities and, 
therefore, should be used. Since the GTD current expressions are 
less complicated than the PO/PTD components and are just as 
accurate, they are sufficient for monostatic work. 

The equivalent currents predictions agree quite well with MM 
and experimental results up to 45° from normal incidence. As 
grazing incidence is approached, the equivalent currents patterns 
deviate from the MM and experimental. These discrepancies can 
probably be resolved by including corner diffraction [4] and 
higher-order terms [14]. 

The method of equivalent currents assumes that each edge is 
an infinite wedge; therefore, it does not properly account for the 
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joining of two edges at a corner. The need for some method of 
including the presence of the corners is illustrated in the 
results for a tilt angle of 0' =45 (Figs. 6 and 7), which 
corresponds to a plane of diffraction through both corners. These 
results deviate from the MM and experimental data to a greater 

o 

extent than for 0' =30 . No rigorously formulated corner 
diffraction coefficients currently exist. Sikta uses 
semi-empirical corner diffraction coefficients [4] and obtains 
accurate results despite the nonuniqueness of the fields in 
certain directions [14]. If carefully applied so that problem 
areas are avoided, this corner diffraction coefficient may lead to 
improved pattern results. 

For vertical polarization higher-order diffractions dominate 
near grazing. The need to include higher-order terms is evident 
in comparing the vertical polarization results (Figs. 5 and 7) to 
the horizontal polarization results (Figs. 4 and 6). The results 
near grazing are much more accurate for horizontal than for 
vertical polarization. Michaeli formulates equivalent currents 
for second-order diffractions [14] that should improve the results 
near grazing incidence. 
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IV. Conclusions 


The scattering patterns that can be accurately calculated 
using GTD, the traditional high-frequency prediction technique, 
are limited to principal planes. This work used the method of 
equivalent currents to extend scattering prediction to 
nonprincipal planes for monostatic RCS calculations. Results 
compared favorably with MM and experimental data away from grazing 
incidence. The addition of corner diffraction and higher-order 


terms was suggested to improve results near corners and grazing. 


V. Publications 


During this reporting period we had the following 

publications : 

a. T. Griesser, C. A. Baianis, K. Liu, "RCS analysis and 
reduction for lossy dihedral corner reflectors, " Proc . IEEE, 
May 1989 (special issue on RCS). 

b. T. Griesser and C. A. Baianis, "Reflections, diffractions, and 
surface waves for an interior impedance wedge of arbitrary 
angle," accepted for publication in IEEE Trans. Antennas 
Propagat . 

c. L. Polka and C. A. Baianis, "Nonprincipal plane scattering 
from perfectly conducting, flat, rectangular plates," 
submitted for presentation at the 1989 International IEEE 
Symposium on Antennas and Propagation. 
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VI . Future Work 


The model for the perfectly conducting flat plate will be 
completed by adding corner diffraction and second-order terms. 
The equivalent currents model will be extended to bistatic 
scattering as well as to lossy plates. Finally, grazing incidence 
scattering will be examined more closely in an attempt to predict 
scattering at and near grazing as a function of tilt angle and 
plate size. Interesting trends have been experimentally observed 
for grazing incidence [3], but no rigorous theoretical models 
exist to match these results. 


- 22 - 


References 


[ 1 1 J. B. Keller, "Geometrical theory of diffraction," J. Opt. 
Soc. Amer., vol . 52, no. 2, pp. 116-130, Feb. 1962. 

[2] R. G. Kouyoumjian and P. H. Pathak, "A uniform geometrical 
theory of diffraction for an edge in a perfectly conducting 
surface," Proc. IEEE, vol. 62, no. 11, pp. 1448-1461, Nov. 
1964. 

[3] R. A. Ross, "Radar cross section of rectangular flat plates 
as a function of aspect angle," IEEE Trans. Antennas 
Propagat., vol. AP-14, no. 3, pp. 329-335, May 1966. 

[4] F. A. Sikta, W. D. Burnside, T. Chu, and L. Peters, Jr., 
"First-order equivalent current and corner diffraction 
scattering from flat-plate structures," IEEE Trans. Antennas 
Propagat., vol. AP-31, no. 4, pp. 584-589, July 1983. 

[5] C. E. Ryan, Jr., and L. Peters, Jr., "Evaluation of 

edge-diffracted fields including equivalent currents for the 
caustic regions," IEEE Trans. Antennas Propagat. , vol. 
AP-17, no. 3, pp. 292-299, May 1969, (also corrections to 

this paper on p. 275, Mar. 1970) 

[6] R. F. Millar, "An approximate theory of the diffraction of 

an electromagnetic wave by an aperture in a plane screen," 

Proc. Inst. Elec. Eng., vol. 103, pt. C, pp. 177-185, Mar. 
1956 (first published as Monograph 152R, Oct. 1955). 

[7] C. E. Ryan, Jr., and R. C. Rudduck, "Radiation Patterns of 
Rectangular Waveguides," IEEE Trans. Antennas Propagat., 
vol. AP-16, no. 4, pp. 488-489, July 1968. 

[8] C. A. Mentzer, L. Peters, Jr., and R. C. Rudduck, "Slope 

diffraction and its application to horns," IEEE Trans. 
Antennas Propagat., vol. AP-23, no. 2, pp. 153-159, Mar. 
1975. 

[9] A. Michaeli, "Equivalent edge currents for arbitrary aspects 
of observation," IEEE Trans. Antennas Propagat., vol. AP-32, 
no. 3, pp. 252-258, Mar. 1984, (also corrections to this 
paper on p. 227, Feb. 1985) 

[10] A. Michaeli, "Elimination of infinities in equivalent edge 
currents, part II: physical optics components," IEEE Trans. 

Antennas Propagat., vol. AP-34 , no. 8, pp. 1034-1037, Aug. 
1986. 


- 23 - 


[11] A. Michael!, "Elimination of infinities in equivalent edge 

currents, part 1: fringe current components," IEEE Trans. 

Antennas Propagat., vol . AP-34, no. 7, pp. 912-918, July 

1986. 

[12] G. T. Ruck, D. E. Barrick, and W. D. Stuart, "Ch. 2 - 

Theory," in Radar Cross Section Handbook, vol. 1, G. T. 
Ruck, Ed. New York: Plenum Press, 1970. 

[13] C. A. Balanis, Antenna Theory: Analysis and Design. New 

York: John Wiley & Sons, Iric., 1982. 

[14] A. Michaeli, "Equivalent currents for second-order 

diffraction by the edges of perfectly conducting polygonal 
surfaces," IEEE Trans. Antennas Propagat., vol. AP-35, no. 
2, pp. 183-190, Feb. 1987. 







o 

o ' 

o 

o 

o 

6 

6 

6 

6 

6 

CM 


CN 

1 

l 

(D 

1 


ujsaa 

ni soa 

OLLVlSONOn 



- 28 - 


90.0 - 45.0 0.0 45.0 90.0 

OBSERVATION ANGLE (THETA) 





-GTD EQUIVALENT CURRENTS EXPERIMENTAL 

-PTD EQUIVALENT CURRENTS 

■MOMENT METHOD 2a=2b=17.18 cm 

f=l 0 GHz 
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-GTD EQUIVALENT CURRENTS EXPERIMENTAL 

-PTD EQUIVALENT CURRENTS 

•MOMENT METHOD 2a=2b=17.18 cm 
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PART B 


PATTERN CONTROL OF HORN ANTENNAS 
I . INTRODUCTION 

Antenna pattern synthesis or the equivalent problem of 
pattern control has been an attractive research topic in the past 
three decades [1 ]-[4] . Particular interest is in the pattern 
synthesis of antenna arraysf 3] , [6] , [7] . The pattern synthesis of 
an aperture antenna has not yet been found to have a satisfactory 
solution . 

Early attempts in using the variational approach to 
optimizing the radiation pattern have been made by Harris et. al . 
[5]. In [5], the optimization of the radiation pattern at a 
single direction is resolved, and a conjugated Green's functional 
distribution for the directivity is suggested. Another attempt in 
optimizing the aperture antenna amplitude pattern was made by 
Elliott [6 ] — [ 8 ] . He discretized the continuous aperture 
distribution, which allowed him to optimize with respect to each 
element. This resulted in a numerical solution to the problem. 

In this report, a functional analysis approach is introduced 
to the optimization, in a range of directions, of the antenna 
amplitude radiation pattern. A continuous eigenfunction solution 
to this problem is found. A demonstration of the use of the 
eigenfunction is also included. 
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II. FAR-ZONE ELECTROMAGNETIC FIELDS 


A general method of finding the radiation pattern of an 
antenna can be obtained by considering the contribution of its 
physical electric current distributions or the equivalent electric 
and magnetic currents. Without loss of generality, we start with 
the magnetic field pattern resulting from the distribution of 
the electric current on a planar structure. 

V,; 

Given the current distribution J(x'.y'), the magnetic field 
can be expressed as: - 


ii.-i-V.A 



1/2 

where k = is the wave constant, and e are the 

permeability and permitivity of the medium, respectively. 

Since we only deal with the far-field pattern, we can assume 
that the observations are made in the far-field of the radiating 
object. Thus we can write that 

R = r - x' sin0 cos <P - y'sinO sin$ (2) 

where r is the distance between the observation point and the 
origin of the coordinate. Since we are only concerned with the 
radiated field component, we can represent the operator V by 

V « > -j k. (3) 
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Using (3) reduces (1) to 


-jkr 


H= 


4nr 


(~jk) x 


, , . jksin0 (x' c.os(p +y' sin <f>) . , . , , „ . 

.J(x ,y )e J J 'dx dy (4) 


We consider a relatively simple case where we observe the 
field along the plane defined by <t> - 0. In addition, we assume 
that the electric current is directed along the X axis. Thus 

k x J = k J(x' ,y' )cos0 a (5) 

h 

where a h is a unit vector aligned with the total magnetic field. 
Finally, combining all the constants in (4) into 1^, the magnetic 
field can be written as 


H[I,0] = a h H 0 cos0 


. jkx' sin0 , 
I(x ) e dx 

*1 


where 


I(x') = 


r y2 

J(x' , y' ) dy' 
Yi 


( 6 ) 


(6a) 


k 

4nr 


(6b). 


The radiated power density pattern is widely used to evaluate 
the quality of an antenna. In the far zone, the radiated power 
density is proportional to the square of the magnitude of the 
radiated magnetic field. Thus we can write that 

P(0) = a | H[ 1 ,0] | 2 (7) 
where a is a proportionality physical constant. 
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III. VARIATIONAL APPROACH TO PATTERN CONTROL 


3. The Construction of the Functional 

The optimization of the radiation pattern in a certain range 
of directions ( from to 0 2 for example), turns out to be a 
variational calculus problem having the following functional, 


JlKx' )] = 


'O'. 


0i 


g(0) P(0) (10 


( 8 ) 


where g(0) is a function which can be used to weight the different 
parts of pattern contributed to the functional. P(0) is 
represented by (7), and it can be expressed as 

P(0) = a | H[I ,0] | 2 = a H [ I ,0 ] • H* [1,0] 


a 


H 0 cos0 


, v jkssin0 , 

1 (s) e ds 

x l 


* 

Ho cos0 

X Y<t) e- Jktslne dt 

V 

X 1 


P(0) = a |h o | 2 cas 2 0 


*2 

*1- 


X ?<s) l*(t) . Jk.lneo-t)^ 

Xl 


(9) 


Substituting (9) into (8) leads to 


J[I(x' )] = 


l 6n 


g(0)«|H o | 2 cos 2 0 


0i 


x 2 

LJxi. 


X 2 


dsdt 


xi 


d0 (10). 


The function inside the middle brackets is a well behaved 
continuous function with respect to 0. Therefore, we interchange 
the order of the two integrations; and for simplicity, let 
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g(0) = l/(a | Ho | 2 ) . Choosing g(0) to be constant implies that all 
parts of the radiation pattern are weighted equally. However, 
g(0) may be chosen as a function of 0 to weight some directions of 
the radiation pattern more heavily in the functional (8), which 
results in a higher optimization level for the given directions. 
After some obvious simplifications, we have: 


J[I(x')] = 


*2 

*1 


*2 * 
1(S)1 (t) 
Xl 


'02 


2 „ jksinO(s-t) .„ 
cos 0 e d0 


0i 


dsdt (11) 


or 




J [I(x')]= 

where 

K(s,t) = 


|X 2 

Xl 


X 2 * 

K ( s , t ) I(s) I (t) dsdt 
Xl 


(12) 


ro 


2 2 n jksin0(s-t) 
cos 0 e 


d6 


I0i 


(12a) 


Notice that K(s,t) is a self-adjoint kernel which satisfies the 
following equation: 


He 


K(s , t) 

= K (t.s) 





(13) 

If the 

functional 

in 

(12) 

should 

have an 

extremal , the 

variational 

derivative 

of 

(12) 

must 

satisfy 

the following 


condition: 

for each point s 0 in the 
source I(s) 

which is equivalent to 


s.Tm 

6 1 s 0 
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0 


(14a) 


111 

h l 


g.Tti] 

6 h 2 



(14b) 


where hj , h 2 are the real and imaginary variational components of 
the current distribution I(x' ), respectively. The variational 
derivative (12a) can be written as 


g.Tti] 

6h t 


1 im 
e-*o 


'x 2 

'x 2 r 


K(s,t) 

J Xp 

JEl L 


I(s) 11 (t) - I (s )1 (t) dsdt 


'x 2 

xi 


( X , s 0 


)I(x)+K(s c ,x)I (x) 


dx= 0 


where I(x) = l(x) + eS(x-s 0 ) 


Following a similar procedure, (12b) gives 


(15). 


K ( s c , x ) I (x )-K(x , s Q ) I (x)j dx = 0 (16) 

Combining (15) and (16), we have an equation which satisfies the 
variational condition of 
'*2 

K(x, s D ) I (x)dx = 0 (17) 

J x 5 

We can define (17) as the extremal condition. 

2. Spectral Analysis of the Extremal Condition 

The solvability of the problem for the extremal condition can 
be analyzed by Fourier transforms. Since K(x,s 0 ) is a function 


AL 


Sh : 


m = 

? So 


X 2 


*1 
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of the difference of the two arguments [ see (12a) ], (17) is 

actually a convolution integral. Thus, the Fourier transform of 
(17) turns out to be 


K( X ) • I ( X ) = 0 


(18) 


where 


I( X ) 


x 2 


1(5) e jX ^ dS 


(19) 


*1 


K( X ) = 


[ X2 jxr 

p 



J x a 

v 


cos 0 e 


-jkrsinO 


d0 


0 1 


dr 


( 20 ) 


To have a non-zero solution to (16), we require that the 
spectrum of 1(X) is perpendicular to the spectrum of K(X). That 
is to say, I (X ) has only a spectral component where K (X ) is zero. 
To find these eigenvalues and their eigenfunctions, we need to 
find the zeros of equation (20). Interchanging the order of 
integration, the inside integral can be evaluated. After some 
simplifications, we have 


K(X) = 


'02 P sin 
cos 0 

01 


(X-ksin0 ) (x^-xt ) 


ksin0 


^ jj ~(*-ksin0) (x 2 +X! )j 


d0 (19) 


Without loss of generality, we pick the origin of the 
coordinate at the center of the aperture. Thus, is zero. 

Also, we denote half of the width as ii). Then (20) can be 
rewritten as 


K ( X ) = 2 


cos 


1 0i 


z a sin{ (X-ksi 00)11)] 

9 X - ksin0 d6 


( 22 ) 
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This can be evaluated numerically. A FORTRAN program was 


written to do the task. Figure 1 shows the spectrum of K(X) when 
plotted with - 5X, = n/16, and G 2 =n/2. Using a numerical 

root finder to determine the zeros of (22), a set of eigenvalues 
can be found. Those are denoted as 

X t , X 2 , , Xi , (23) 

and their corresponding eigenfunctions are 

e" jXlX ', cf jX2X '.---, e- jXiX \ (24) 

The current distribution on the aperture, which satisfies the 
extremal condition can be expressed as linear combinations of 
these eigenfunctions. 
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IV. CONCLUSIONS 


1. Solutions 

From the work of III, we conclude that the extremals of the 
control functional can be written as 

I(x' )= J Ci e" jXiX ' (25) 

1 = 1 

and the corresponding radiation pattern can be expressed as 
P(0) « a 

' i 

where \i * s are the eigenvalues of (22) in 111; C L ' s are the 
spectra of each component. We found that (25) and (26) take a 
surprisingly simple form in their expressions. Although in the 
case when the source distribution is of the finite extent, the 
eigenfunctions, obtained by such a procedure, still play an 
important part. The Ritz method [9] can be used to refine the 
extremal, and to optimize the Cj ' s used in (25) and (26). 

2. An Example 

To demonstrate the derived algorithm, a 5\ aperture antenna 
was examined. The objective was to optimize the antenna pattern 
to allow the largest possible percentages of radiated energy into 

o o o 

the main beam ( 0 ^ 0 :£ 11.25 where 0 = 0 is normal to the 

o o 

apperture ) while reducing the minor lobe ( 11.25 £ 0 £ 90 ) 

energy to the smallest level. Control directivities are chosen to 
be from 71/16 to 7t/2 ( since this is only a half space antenna 

problem ) . 


0 sinKXj -ksin0)u ] 
Ll COS ° Xj - ksinQ 


(26) 
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Figure 1 displays the spectrum of K(X) for this particular 
design. Notice that the zero, which is the closest to X=0, has an 
eigenvalue of 0.4587. Figure 2 shows three normalized (to their 
maximum values) radiation patterns of the eigenfunctions which 
correspond to the first three smallest eigenvalues (in magnitude). 

Figure 3 exhibits a comparison of the two 
pattern-optimization approaches. The solid line is the pattern 
obtained by using the Green's functional distribution suggested by 
Harris et. al[5]; it maximized the radiation pattern only at 0=0. 
The dashed line is the optimization pattern obtained by the method 
that was outlined here. It reduced the percentages of the 
radiated energy from tc/ 16 to n/2 to a lowest possible limit. 
Doing this, the beam efficiency of such an antenna is greatly 
enhanced. Moreover, such a radiation pattern shows a certain 
flatness along the main radiation beam which most antenna 
designers would like to have. The electric current distribution 
chosen for this design is the one which has the smallest 
eigenvalue (the first eigenpattern in Figure 2). 

Figure 4 displays the validity of the method. The solid-line 
shows the ratio variation, with respect to the eigenvalue, of the 
energy radiated within the controlled directions to the main-beam 
while the dashed-line shows the variation of the normalized 
radiated power with respect to the eigenvalue. The curve shows a 
minimum ratio at the eigenvalue that we have chosen. Also, from 
the comparison of the two curves, our algorithm does choose the 



best design from the trade-off of the beam efficiency and the 
maximum radiation power from the main beam. 
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V. FUTURE WORK 


4 


Future work on this project, will concentrate on the 
fo] lowing: 

1. The formulation and the solution of the two-dimensional 
(or three-dimensional, if needed) functional source 
distribution for horn antenna pattern control. 

2. Studies of resistive tapering and its effect on the 
aperture field distribution of the horn antenna, and on 
the overall far-zone radiation pattern. 

3. Numerical solution of the exact boundary value problem of 
the horn antenna. This step is also very important 
because we need it to demonstrate the resultant optimum 
pattern . 
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Figure 3. Radiation pattern of Green’s functional distribution 
and the optimum distribution Control between 1 1 .25 to 90 degrees 
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Figure 4. Radiated power of Main beam and minor lobes 



